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Abstract 

By use of the two measures presented recently, the indivisibility and the backflow of informa- 
tion, we study the non-Markovianity of the dynamics for a two-level system interacting with a 
zero-temperature structured environment without using rotating wave approximation (RWA). In 
the limit of weak coupling between the system and the reservoir, and by expanding the time- 
convolutionless (TCL) generator to the forth order with respect to the coupling strength, the 
time-local non-Markovian master equation for the reduced state of the system is derived. Under 
the secular approximation, the exact analytic solution is obtained and the sufficient and necessary 
conditions for the indivisibility and the backflow of information for the system dynamics are pre- 
sented. In the more general case, we investigate numerically the properties of the two measures for 
the case of Lorentzian reservoir. Our results show the importance of the counter-rotating terms to 
the short-time-scale non-Markovian behavior of the system dynamics, further expose the relations 
between the two measures and their rationality as non-Markovian measures. Finally, the complete 
positivity of the dynamics of the considered system is discussed. 
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I. INTRODUCTION 



Realistic quantum systems cannot avoid interactions with their environments, thus the 
study of open quantum systems is very important. It is not only relevant for better un- 
derstanding of quantum theory, but also fundamental for various modern applications of 
quantum mechanics, especially for quantum communication, cryptography and computa- 
tion The early study of dynamics of open quantum systems usually consists in the 
application of an appropriate Born-Markov approximation, that is, neglects all the memory 
effects, leading to a master equation which can be cast in the so-called Lindblad form 0, [J. 
Master equations in Lindblad form can be characterized by the fact that the dynamics of 
the system satisfies both the semigroup property and the complete positivity, thus ensuring 
the preservation of positivity of the density matrix during the time evolution. We usu- 
ally attribute the dynamical processes with these evolutional properties to the well-known 
Markovian ones. 

However, people recently found that Many relevant physical systems, such as the quan- 
tum optical system [4|, quantum dot (5)], superconductor system 6], could not be described 



simply by Markovian dynamics. Similarly, quantum chemistry W and the excitation trans- 
fer of a biological system jg| also need to be treated as non-Markovian processes. Quantum 
non- Markovian processes can lead to distinctly different effects on decoherence and disen- 
tanglement (9), [lo| of open systems compared to Markovian processes. These non-markovian 
effects can on the one hand enrich the basic theory of quantum mechanics, on the other 
hand benefit the quantum information processing. Because of these distinctive proper- 
ties and extensive applications, more and more attention and interest have been devoted 
to the study of non-Markovian processes of open systems, including the measures of non 
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of non-Markovian environment has been realized. 

The measure of non-Markovianity of quantum evolution is a fundamental problem which 
aims to detect whether a quantum process is non-Markovian and how much degrees it de- 
viates from a Markovian one. Based on the distinguishability of quantum states, Breuer, 
Laine and Piilo (BLP) 11] proposed a measure to detect the non-Markovianity of quan- 
tum processes which is linked to the flow of information between the system and environ- 



ment. Alternatively, Rivas, Huelga and Plenio [13J (RHP) also presented a measure of 
non-Markovianity by employing the dynamical divisibility of a trace-preserving completely 
positive map. It is clear that the BLP measure is based on the physical features of the 
system-reservoir interactions, while the RHP definition is based on the mathematical prop- 
erty of the dynamical maps. It has been shown that the two measures agree for several 



important and commonly-used models [33J, but do not agree in general 34]. In this paper, 
we will use both the two measures to describe the non-Markovianity of the dynamics of the 
considered system, so as to more clearly see their relation, as well as the rationality as the 
measure of non-Markovianity. 

The study of the dynamics of non-Markovian open quantum systems is typically very in- 
volved and often requires some approximations. Almost all the previous treatments are based 
on the RWA, that is, neglect the counter-rotating terms in the microscopic system-reservoir 
interaction Hamiltonian. However, the counter-rotating terms which are responsible for the 
virtual exchanges of energy between the system and the environment not always can be 
neglected. For example, for the wide-frequency-spectra reservoir or when the frequency dis- 
tribution of the structured environment is detuned large enough from the transition of the 
system, the RWA is invalid. Another motivation of this paper is thus to study the effect of 
the counter-rotating terms on the non-Markovian dynamics of the considered open quantum 
system. 

The article is organized as follows. In Sec. II we introduce the microscopic Hamiltonian 
model between the system and its environment, and derive the non-Markovian time- local 
master equation for a two-level system weakly coupled to a vacuum reservoir, by using the 
TCL approach to the forth-order but without employing RWA in the interaction Hamilto- 
nian. In Sec. Ill, we investigate the non-Markovianity of the system dynamics in terms of 
both the RHP and BLP measures. Through the analytical solution in the secular approx- 
imation, we obtain the sufficient and necessary conditions for the dynamical indivisibility 
and the backflow of information, showing the effect of the counter-rotating terms on the 
non-Markovian dynamics of the system, and exposing the relations between the BLP and 
RHP measures. In sec. IV, by choosing the Lorentzian spectra reservoir as an exemplary 
example, we further demonstrate the effect of the counter-rotating terms on the dynamical 
indivisibility and the backflow of information, and clarify the rationality of the two non- 
Markovian measures. Finally in Sec.V, we discuss simply the complete positivity of the 



system dynamics. And the conclusion is arranged in Sec. VI. 



II. THE MICROSCOPIC MODEL 

Consider a two-level atom with Bohr frequency uq interacting with a zero-temperature 
bosonic reservoir modeled by an infinite chain of quantum harmonic oscillators. The total 
Hamiltonian for this system in the Schrodinger picture is given by 

H = -u Q a z + ^u k blb k + ^# fc (cr + + o_)(b k + b%), (1) 
1 k k 

where a z and a± are the Pauli and inversion operators of the atom, u k , b k and b^ are respec- 
tively the frequency, annihilation and creation operators for the k-th harmonic oscillator of 
the reservoir. The coupling strength g k is assumed to be real for simplicity. The distinct 
feature of this Hamiltonian is the reservation of the counter- rotating terms, cr + b^ and <J~b k , 
which is the so-called without RWA we call in this paper. Note however that our starting 
point is the dipole interaction Hamiltonian between the atom and its environment, whose 
derivation starting with the canonical Hamiltonian involves the discarding of a term which 
is quadratic with respect to the radiation field. The discarding is not based on the RWA, 
but the fact that for low-intensity radiation, the quadratic term is much small compared to 



the dipole interaction one 36]. 



The time-convolutionless projection operator technique is most effective in dealing with 
the dynamics of open quantum systems. In the limit of weak coupling between the system 
and the environment, by expanding the TCL generator to the forth order with respect 
to coupling strength, the non-Markovian master equation describing the evolution of the 
reduced system, in the interaction picture, can be written as [For the main clue of its 
derivation, see appendix A.] 

^ = -i[H LS (t), p(t)\ + D[p(t)\ + D'[p(t)\, (2) 

where 

H LS {t) = S + (t)a + a^ + S-{t)<r-a+, (3) 
is the Lamb shift Hamiltonian which describes a small shift in the energy of the eigenstates 



of the two-level atom. In many theoretical researches 



231 ]. this term was neglected usually. 



But in this paper, we will take it into the consideration. The Lamb shift includes the second 
and forth order contributions, 

S ± (t) = Si I (t) + S?(t), (4) 

which respectively come from the second and forth order perturbative expansion of the TCL 
generator. The second order Lamb shift is 

S£(t) = ±J*dr J duoJ{uj) sin[(w =F w)r], (5) 

with J{oj) = J2k \9k\ 2 fi(w — ^k) the spectral distribution of the environment. The expression 
for the forth order Lamb shift S ! ± (t) is cumbersome which is presented in the appendix A. 
The dissipator D[p(t)] that describes the secular motion of the system has the form 

D[p(t)} = T_(t){a_p(t)a + - ±{a + a_,p(t)}} (6) 

+ T + (t){a + p(t)a_ - ±{a„a + ,p(t)}} 

+ T (t){a + a^p(t)a + a^ - p(t)}}, 

where the first line describes the dissipation of the atom to the vacuum environment with 
time-dependent decay rate and the second line denotes the heating of the atom in the 

vacuum environment with time-dependent heating rate This heating is related to the 

dissipation, for a ground-state atom in a zero-temperature environment, there is no heating 
effect. Dissipation and heating are usually accompanied by decoherence. The last line in 
eq.(6) describes the pure decoherence with time-dependent decoherence rate T (t). The 
time-dependent transition rates r±(i) also include the second and forth order perturbative 
contributions of the TCL generator, 

r ± (0 = r^) + rm (7) 

with the second order contribution as 

r£ 7 (t) =2^'rfr J du>J(u)co8[(u>o±w)T\. (8) 

While r (t) completely comes from the forth-order perturbative contribution. All the forth- 
order contributions are presented in the appendix A. Eq.(6) indicates that the dissipative 
model of eq.(l), except for inducing the energy exchange between the system and its envi- 
ronment, also makes decoherence of the system. But the rate of decoherence is much less 
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than that of energy dissipation, because To(t) is only a forth-order contribution term of TCL 
perturbative expansion. 

The dissipator D'[p(t)} represents the contribution of the so-called nonsecular terms, that 
is, terms oscillating rapidly with Bohr frequency Uq, 

D'[p(t)] = [a(t) + iP(t))a + p(t)a+ + h.c, (9) 

here h.c. denotes the Hermitian conjugation. These nonsecular terms sometimes may also be 
neglected under the so-called secular approximation [35|. The time-dependent coefficients 
a(t) and (3{t) also include the second and forth order contributions, 

a(t) =a H (t) + a IV (t), (10) 

P(t)=(3 n (t)+(3 IV (t), (11) 

with 

a n (t) = 2 J* dr J duJ(u) cos[u(t - r)} cos[u {t + r)], (12) 

and 

cos[u{t — r)] sin[o;o(t + r)]. (13) 

The forth-order contributions are listed in the appendix A. 

Note that the dynamics for a two-level system embedded in a zero-temperature structured 
environment, under RWA, can be solved exactly, where the corresponding master equation 



has the Lindblad-like form 



j/(t) = - l -S{t)[a + a^p{t)]+ 1 {t){a^ P {t)a + - \w + °-, p(t)}}, (14) 

where the time-dependent decay rate 7(t) and Lamb shift S(t) are related to the correlation 
function of the reservoir. Comparing this equation with eq.(2), we see that the last two 
terms in the dissipator D[p(t)], that is, the heating and the pure decoherence terms, as 
well as the nonsecular dissipator D'[p(t)} and the Lamb shift 5*_(t), are completely from the 
contribution of the counter-rotating terms presented in the interaction Hamiltonian. While 
the decay rate T_(t) and the Lamb shift S+(t) include the contributions of both rotating 
and counter- rotating terms, but the main contributions [i.e., the second-order terms r{. 7 (t) 
and St (t)] come from the rotating terms. In fact, by expanding the decay rate ^(t) and 
the Lamb shift S(t) to the second order with respect to coupling strength, one obtain r{_ 7 (t) 



1, 



and S+(t) |4(. In the following, we will show that the contributions that come from the 
counter-rotating terms are important, in particular to the short-time-scale non-Markovian 
behaviors. 

III. MEASURES OF NON-MARKOVIANITY 

Recently, people have been interested in the study of non-Markovianity of open quantum 
systems. Several definitions or measures of non-Markovian dynamics have been 



presented. In this section, we will employ two of the measures, i.e., the RHP [3] and BLP 



ll| measures, to investigate the non-Markovian dynamics of the considered system so as to 
see the effect of the counter-rotating terms on non-Markovianity and the relation between 
the two measures. 

A. Divisible and indivisible dynamics 

A trace-preserving completely positive map e(t 2 ,0) that describes the evolution from 
times zero to t 2 is divisible if it satisfies composition law, 

e(t 2 ,0)=e(t 2 ,ti)e(ti,0), (15) 

with e(t 2 ,ti) being completely positive for any t 2 >t\> 0. Due to the continuity of time, 
eq.(15) is always fulfilled in form. The key point for divisibility is actually the complete 
positivity of e{t 2l ti) for any t 2 > t± > 0. If there exist times t\ and t 2 such that the map 



e(t 2 ,t±) is not completely positive, then the dynamical map e(t 2 ,0) is indivisible. RHP (13 1 
defined all the divisible maps to be Markovian. Therefore, the indivisibility of a map ad- 
vocates its dynamical non-Markovianity. It was shown that all the evolutions governed by 
Lindblad-type master equation with positive transition rates are divisible jszj], thus Marko- 
vian. 



It was proved 13j that the indivisibility of map e(t, 0) is equivalent to the complete 



positivity of the quantity, 

(t) = lim \\Ht + e j)® mm -^ 

Only for divisible map, g(t) = 0. Where |$) is a maximally entangled state between the 
system of interest and an ancillary particle, and the map e performs only on the state 



of the system. Using the time- local master equation ^ = C t (p), this expression may be 
equivalently written as k| 

= ||[/ + (£,®/)c]|$)($|||-l 

v 1 e^o+ e v ' 

The function g(t) is the so-called RHP non-Markovian measure. If and only if g(t) = 
for every time t e {O,^}, the map e(t2,0) is Markovian. Otherwise it is non-Markovian. 
The distinctive advantage of RHP non-Markovian measure is that its calculation can be 
processed only by the use of time-local master equation, not requiring the exact form of the 
dynamical map e(t, 0). In the following, we call the time interval that satisfies g(t) > the 
indivisible dynamical interval (IDI). For a non-Markovian process, there must exist one or 
several IDIs. 

For the open two-level system considered in this paper, suppose that |$) = ^[|01) + |10)], 
a straightforward deduction using equations (2) and (17) gives 



g = -|r_ + r + + ^(r_ - r + ) 2 + 4 (a 2 + /? 2 )| (is) 
+ ||r_ + r + - ^(r_-r + ) 2 + 4(a2 + /3 2)| 

+ i[|r |-r -2r_-2r+], 

where for compactness we omit the argument of all the time-dependent coefficients. Obvi- 
ously, the Lamb shift Higit) has no effect on the indivisibility of the system dynamics. 



B. Backflow of information 



The second measure of non-Markovianity for quantum processes of open systems we 



employ is proposed by BLP [llj] which is based on the consideration in purely physics. Note 
that Markovian processes always tend to continuously reduce the trace distance between 
any two states of a quantum system, thus an increase of the trace distance during any 
time interval implies the emergence of non-Markovianity. BLP further linked the change 
of the trace distance to the flow of information between the system and its environment, 
and concluded that the back flow of information from environment to the system is the key 
feature of a non-Markovian dynamics. In quantum information science, the trace distance 
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for quantum states p\ and P2 is defined as [l| 

D(pi,p 2 ) = ~tr\pi - p 2 \, (19) 

with \A\ = V A + A. For a given pair of initial states pi,2(0) of the system, the change of the 
dynamical trace-distance can be described by its time derivative 

a(t,p lt2 (0)) = j t D( Pl (t),p 2 (t)), (20) 

where pi^{t) are the dynamical states of the system with the initial states /^(O). For Marko- 
vian processes, the monotonically reduction of the trace distance implies cx(i, pi^O)) < 
for any initial states pi,2(0) and at any time t. If there exists a pair of initial states of the 
system such that for some evolutional time t, a(t, pi^O)) > 0, then the information takes 
backflow from environment to the system, and the process is non-Markovian. 

In order to calculate the BLP measure, we must solve the dynamics of the system. For 
this purpose, we write the alternative Bloch equation of eq.(2) as [see appendix B for their 
derivation], 



K = --(r_ + r + + r - 2a)b x + (sl -s+- f3)b y , (21) 

k = -i(r_ + r + + r + ia)b y - (sl - s + + p)b x , (22) 
& 2 = -(r_ + r + )6 2 + r + -r_, (23) 

where the three components of the Bloch vector are defined as bj(t) = Tr[p(t)<jj} with 
j — x,y,z and <7j the Pauli operators. In terms of Bloch vector, the trace distance of eq.(19) 
may be expressed as 

D(t) = \sj{Ab x f + {Ab y f + (Ab z y (24) 

where Abj = bij(t) — b 2 j(t) are the differences between the two Bloch components at evolu- 
tional time t. Correspondingly, the derivative of this trace distance becomes 

a = -\[{Ab x ) 2 + (Ab y ) 2 + (A6 2 ) 2 ]- 1/2 {(r- + T + + r - 2a)(Ab x ) 2 (25) 
+ (r_ + T + + r + 2a){Ab y ) 2 + A/3(Ab x )(Ab y ) + 2(r_ + r + )(A6,) 2 }, 

where we have used the Bloch eqs.(21)-(23) in the deduction process. According BLP's 
criterion, o > indicates the backflow of information from environment to the system. In 
the following, we call the time intervals in which a(t) > the information-backflow intervals 
(IBIs). 

9 



C. Secular approximation 



In order to see the effect of counter-rotating terms and make a distinct comparison be- 
tween the BLP and RHP measures in the current system, we now consider the case where 
the nonsecular term D'[p(t)} can be neglected, i.e., performing the so-called secular approxi- 
mation. Here for the sake of discrimination, we call as in many literatures the rotating-wave 
approximation that used after tracing over the bath degrees of freedom the secular approxima- 
tion. In other words, the secular approximation and the RWA have the same mathematical 
approaches-throwing away the rapidly oscillating terms in time, merely the times the ap- 
proximations taking place are different. Just as pointed out in reference {35 ], this kind 
of secular approximation though also is an average over rapidly oscillating terms, it does 
not wash out the effect of the counter-rotating terms present in the coupling Hamiltonian. 
Under the secular approximation, the master equation (2) has the Lindblad-like form with 
time-dependent transition rates, r±(i), Vo(t) and Lamb shift HLs(t). Employing the method 
proposed in [39}, the corresponding Bloch eqs.(21)-(23) in this case can be solved exactly 



which gives 



with 



b x (t) = e- (%.(O) cos 5(t) - b y (0) sin 5(t)}, (26) 

b y {t) = e- 0it) [b x {O) sm5{t) + b y (0) cos<?(t)], (27) 

K{t) = e" A « {6,(0) + fjse^[Y + {s) - r_(a)]} , (28) 

©(*) = \f ds[T_{s) + r+(s) + r (s)], (29) 
2 jo 

A(t)= f ds[Y„{s)+T+{s% (30) 
J 



and 



6(t) = I ds\S + {s) - SJs)]. (31) 
J 

Inserting these solutions into eq.(25), we get 

a = -i/(t){e- 20 W(r_ + T + + T )[(Ab x (0)) 2 + (^b y (0)) 2 ] (32) 
+ 2e- 2A W(r_ + r + )(A6 2 (0)) 2 }, 

where I(t) = {e' 20 ^ [(Ab x {0)) 2 + (A^(0)) 2 ] + e- 2A V (Ab^O)) 2 }- 1 / 2 is a positive function 
and Abj(0) = feij(0) — 6 2j (0) is the difference between the two initial Bloch components. This 
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expression shows that the sufficient and necessary conditions for the backflow of information 
from environment to the system are 

r_(t) + r + (t) + r (t) <o, (33) 

or 

r_(t) + r+(t) <o. (34) 

Because if at some time t, one of this conditions is satisfied, then we can always find a pair 
of initial states such that a(t) > 0. For example, if eq.(33) fulfils, it suffices to choose the 
initial states satisfying Ab z (0) = 0. Conversely, if a{t) > at some time t, then at least one 
of the two conditions must be satisfied. 

On the other hand, under secular approximation, eq.(18) is simplified as 

g = i{2|r_| + 2|r + | + |r | - 2r_ - 2r + - r }. (35) 

Obviously, when one of the three rate functions, T_(t), T + (t) or r (t), is negative, then 
g > 0, vice versa. Thus the sufficient and necessary conditions for the indivisibility of the 
dynamics are 

r_(t) < 0,or r+(t) < 0,or r (t) < O. (36) 

Eqs. (33), (34) and (36) demonstrate two important results. One the one hand, the 
counter-rotating terms [which induce r+(i), To(t) and a part of r{_ y (t)] may have important 
effect to the non-Markovian dynamics of the system, according to RHP and BLP measures. 
On the other hand, they show that the conditions for the backflow of information are much 
more rigorous than that of indivisibility. The later only requires one of the transition rates 
to be negative, while the former further requires the sum of the two or the total transition 
rates to be negative. This conditionality once again validates the previous results: The 
backflow of information must lead the indivisibility of the dynamics, but the reverse is not 
true 34]. However, for Lindblad-like master equation with only single transition rate, the 
sufficient and necessary conditions for the two measures become clearly the same, denoting 



the consistency of the two measures in this case [33j. 



IV. NON-MARKOVIAN DYNAMICS FOR LORENTZIAN SPECTRUM 



In order to further demonstrate quantitatively the effect of the counter-rotating terms, 
as well as the rationality of the two non-Markovian measures, we specify our study to a 
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particular reservoir spectra, Lorentzian spectra, 



JM = 2n[(uo-u-Ay + Xr (3?) 



which describes the interaction of an atom with an imperfect cavity and is widely used in 
literatures. Where uq denotes the transition frequency of the atom, A = ojq — u c is the 
frequency detuning between the atom and the cavity mode. A is the width of Lorentzian 
distribution, which is connected to the reservoir correlation time r# = A" 1 . The parameter 
7o can be regarded as the decay rate for the excited atom in the Markovian limit of flat 
spectrum which is related to the relaxation time t$ = 'Jq 1 . For the Lorentzian spectra, all 
the time-dependent coefficients including S±(t), T±(t), T (t), a(t) and f3{t) can be calculated 
analytically, but the expressions are too complicated. We thus study them only numerically. 

In Fig.l, we show the time evolution of these coefficients. For our purpose, we intention- 
ally choose three special sets of parameters. It shows that for narrow spectrum and small 
detuning [In Fig.l (a) and (d), = 0.2%, A/uo = 2%], plays the dominant role, 

while r + (t) and T (t) are almost zero. The nonsecular coefficients a(t) and (3(t) in this 
case behave fast oscillations [Fig.l (d)], so that on average in time the effect can also be 
neglected. These results imply that for this set of parameters, the counter-rotating terms 
in Hamiltonian (1) play little effect actually to the system dynamics and the commonly- 
used RWA is valid. However, for wider spectrum or/and larger detuning, the results are 



different [see Fig.l (b) and (c)], where though T is still near zero [38j], r + clearly can not 
be neglected. Thus the counter-rotating terms in these cases are important and the RWA 
is invalid. Of course, for very wide spectrum, one may expect that the dynamics tends to 
be Markovian. The positivity of the T±(t) and To(t) in Fig.l (c) confirms this point. In 
addition, when A is small, the correlation time of the environment is longer, thus T_ in Fig.l 

(a) oscillates to emerge negative values in a relatively longer time. With the increasing of A, 
the correlation time becomes small and small, and the times for to be negative shorten 
or even vanish [Fig.l (b) and (c)]. Note that the observable negative values of r + in Fig.l 

(b) demonstrate the contribution of the counter-rotating terms to the non-Markovianity of 
the system dynamics. 

Note that in the RWA, the corresponding master equation (14) may be solved exactly. 
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For the Lorentzian spectrum, the RHP and BLP measures may be expressed as 33 j 



and 



where 



7(t) = Re 



m - < for Mt) £ (38) 

-7(t) for 7 (t) < 



a(t) = -j(t)F(t). (39) 
27 Asinh(c/t/2) 



dcosh(dt/2) + (A - iA) sinh(rft/2) 



(40) 



with d = y (A — zA) 2 — 270A. The positive real function F(t) is defined as, 



2 -|r(t) , 1, 12 -ir(t) 

F(t) = - + (41) 
^g-rw + | 6 | 2 

with r(t) = fidt-ytf), and a = <l| Pl (0)|l> - <l|p 2 (0)|l>, 6 = <l| Pl (0)|0> - <l|p 2 (0)|0> being 
the differences of the population and of coherence respectively for the two given initial states. 
Eqs.(38) and (39) show that under the RWA, the distributions of IDIs and IBIs are exactly 
the same, which are determined by j(t) < 0. In the following, we study numerically the 
evolution of the measures a and g, under the condition without using RWA, so as to further 
highlight the non-Markovian effect of the counter-rotating terms, as well as the rationality 
of BLP and RHP measures. 

In Fig. 2, we show the time evolution of the measure a in the same parameters as in 
Fig.l, where the solid lines are plotted according to eq.(25) and the dot lines according to 
eq.(39). We choose the pair of initial states to be pi(0) = |1)(1| and p2(0) = |0)(0|, which 



can maximize the BLP measure III]. For evidence, we only give the time intervals of a > 0, 
i.e., the IBIs. We can see clearly the corrections of the counter- rotating terms on the BLP 
measure. In Fig. 2 (a), both the distributions of the IBIs and the shapes of the two curves 
are similar, responding that the counter-rotating terms make lesser effect to the backflow of 
information in this case which is in line with the idea of RWA. The dips on each peaks of the 
solid-line in Fig. 2 (b) are due to the negativity of r+(t) at that times [see Fig.l (b)], implying 
that T + has the offset on the backflow of information. There is no IBI in Fig. 2 (c), denoting 
that under the choice of this set of parameters, there is no backflow of information, or 
equivalently the dynamics is Markovian according to BLP measure, which is in line with the 
non- negativity of T±(t) and r (t). In addition, the time scale for the backflow of information 
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is consistent with the reservoir correlation time A -1 [Fig. 2 (a), (b)]. All these results show 
that on the one hand the counter-rotating terms can affect the backffow of information, and 
on the other hand the correction of the counter-rotating terms on the backffow of information 
is reasonable. 

In Fig. 3, we plot the time evolution of the measure g in the same parameters as in 
Fig.l. We see that when the counter-rotating terms are omitted, the distribution of the 
IDIs agrees with that of IBIs [see the dot lines in Figs. 2 and 3]. The non-Markovian time 
scale predicted by measure g is also in accordance with the reservoir correlation time A -1 . 
The horizontal dot line in Fig. 3(c) denotes that under the choice of those parameters, the 
dynamics is actually Markovian. All these results show that with no counter-rotating terms, 
the RHP and BLP measures agree. Both of them can depict rightly the non-Markovianity 
of the underlying dynamics. However, when the counter-rotating terms are considered, 
the case is distinctly different: The IDIs now become (0, oo) [see the solid lines in Fig.3], 
which are clearly inconsistent with the practice. Because first of all, the non-Markovian 
time scales in the underlying conditions are never infinite. Next, in Fig. 3(a), the choice 
of the parameters is consistent with the RWA, the result after considering counter-rotating 
terms should has some tiny, not distinct amendments, over the result under RWA. For 
the parameters in Fig. 3(c), the reservoir correlation time A -1 is very short and the system 
dynamics is actually Markovian, should not appearing long-time non-Markovianity. These 
egregious results denote that the RHP measure in these cases is invalid. Note that the reason 
for resulting in these unpractical phenomena is mainly due to the nonsecular coefficients a(t) 
and (3(t). When these nonsecular coefficients are neglected, eq.(35) is not seen to deviate 
obviously from the practice. 



V. COMPLETE POSITIVITY 

The evolution of a real physical state should be not only positive but also complete 
positive. In practical theoretical study, however, due to the application of some assumptions 
and approximations, the positivity or the complete positivity may not always be satisfied. 
Here we present a study of the complete positivity for our considered model, i.e., the master 
equation of (2). As the damping matrix has the block diagonal form (see appendix B), 



thus we can directly use the conditions for complete positivity presented by Hall 



39|. The 
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necessary condition of the complete positivity, for the master equation (2), may be given by 
two inequalities: 

A(t) > 0, (42) 

26(t) > A(t), (43) 

with Q(t) and A(t) given by eqs.(29)-(30). The sufficient condition is also given by two 
inequalities. The first one coincides with eq.(42) and the second one may be expressed as 

X (t) cosh9(t) < 1 + A 2 (t) - K 2 (t) - 2\A{t) - X (t)\, (44) 

where X if) = e~ 2@ ^ , A(t) = e~ A W, n{t) = A(t) fids\T+(s) - T^(s)]A-\s), and 9{t) = 
2 Jq ds^Ja 2 (s) + (3 2 (s) with a(t), f3(t) given by eqs.(10)-(ll). Using inequality (43) to release 
the modulus in the right-hand side, we get 

X (t) cosh9(t) < [1 - A(t)} 2 + 2 X (t) - n 2 (t), (45) 

Note that the left-hand side of eq.(45) is relevant to the nonsecular motion, but the right- 
hand side only depends on the secular motion. As 9(t) increases with time t, eq.(45) is 
not satisfied for long times. But in short non-Markovian time scales we are interested in, 
it may be fulfilled. In order to see this, we plot the time evolution of function G(t) = 
[1 — A(t)} 2 + 2x(t) — K, 2 (t) — x(t) cosh. 9 (t) as in Fig. 4, for the same parameters as in Fig.l 
and under the Lorentzian spectra. Obviously, in the scale of the correlation times A -1 , 
G(t) > 0. The condition of eqs.(42)-(43) is satisfied for all times in this case. Thus in the 
short non-Markovian time scales, the evolution of the system is physical. In the secular 
regime, the sufficient condition eq.(45) can be relaxed to 

[l-A(t)] 2 + X (t)-n 2 (t)>0, (46) 

which can be satisfied for much more longer times for the Lorentzian reservoir. 

VI. CONCLUSION 

In conclusion, we have studied the non-Markovianity of the dynamics for a two-level 
system interacting with a zero-temperature structured environment without using RWA. In 
the limit of weak coupling between the system and its reservoir, by expanding the TCL 
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generator to the forth order with respect to the coupling strength, we have derived the time- 
local non-Markovian master equation for the reduced state of the system. Under the secular 
approximation, the TCL master equation has the Lindblad-like form with time-dependent 
transition rates. We have obtained the exact analytic solution. The sufficient and necessary 
conditions for the indivisibility and the backflow of information for the system dynamics 
were presented, which showed two important results: First, the counter-rotating terms may 
play important roles to the indivisibility and the backflow of information for the system 
dynamics. Second, it showed explicitly that the BLP and RHP measures generally do not 
coincide. It demonstrated more clearly the previous result: The backflow of information 
must lead to the indivisibility of dynamics, but the reserve is not true. 

When the nonsecular terms are included, we have investigated numerically the non- 
Markovian properties of the system dynamics by assuming that the environment spectrum 
is Lorentzian. By compared with the result under RWA, we found that the BLP measure 
is corrected appropriately, but the RHP measure is inconsistent with practice, showing that 
the RHP measure has finite applicable range. 

Finally, we have discussed the complete positivity of the underlying dynamics. We have 
presented the sufficient and necessary conditions of the complete positivity. Numerical 
simulation showed that these conditions can be satisfied in the short non-Markovian time 
scale. 

The measure of non-Markovianity is a fundamental problem in the study of open quan- 
tum system dynamics. Although several measures of non-Markovianity have been presented 
already, it is noted that these measures are not completely equivalent to each other. There- 
fore, the problem for measuring the non-Markovianity of quantum processes still remains 
elusive and, in some sense, controversial. At present stage, it is meaningful and necessary to 
expose the characteristics of various measures and their relations in some concrete systems. 

The investigation of a two-level system interacting with a bath of harmonic oscillators, i.e., 
the spin-boson model, is of particular interest in the theory of open quantum system. In the 
context of quantum computation, it represents a qubit coupled to an environment, which 
can produce dissipation and decoherence. Though in the numerical simulations we have 
only considered the Lorentzian environment, our analytic results adapt to other structured 



environments, such as the Ohmic reservoir, the photonic band-gap material 4Gj, etc. By 



properly engineering the structure of the environment, one can control the non-Markovian 
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dynamics of the open quantum system, so as to effectively control the evolution of some 
interesting physical quantities, such as the quantum coherence, quantum entanglement and 
discord, etc. Therefore, our work will be helpful for the quantum information processing. 

Of course, our model is not fully general. First of all, we have considered only a two- 
level system weakly coupled a zero-temperature environment. Next, our starting point is 
based on the dipole interaction Hamiltonian between the atom and its environment, not 
on the canonical Hamiltonian. Finally, we have used the TCL perturbation expansion for 
the derivation of master equation eq.(2). Thus our results are still conditional and further 
investigations may be necessary. 
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Appendix A: Derivation of the master equation and the time-dependent forth-order 
coefficients 

In our study, the derivation of the forth-order TCL master equation (2) is very cum- 
bersome. Here we can present only the main clue about the deduction. Our calculation is 



based on the description of reference (4] about the TCL projection operator technique. By 
assuming a factoring initial condition p(0) = ps(0) <S> Pb for the system and environment, 



one obtains a homogeneous TCL master equation [see (9.33) of J4]] 

^Vp(t) = K{t)Vp{t). (Al) 

Due to the assumption of vacuum reference state ps = |0)(0| for the environment, the TCL 
generator K.(t) only has even-order terms in its perturbation expansion. The second- and 
forth-order TCL generators may be calculated directly via eqs.(9.61)-(9.62) of reference J4], 
where the related operators and in the interaction picture are given by, 

F k (t) = a + e iulot + a.e~ iul0 \ (A2) 
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Qk(t) = g k (b k e~^ t + bte^ t ). 



(A3) 



Calculating the second- and forth-order TCL generators and sorting them in operators, then 
eq.(Al) reduces to the required master equation. 

In the master equation (2), each of the time-dependent coefficients consists of in principle 
two parts-the second and the forth order parts. The second-order parts have relatively 
simple expressions, but the expressions of the forth-order parts are very complex. In terms 
of abbreviation ty —U — tj with t = t, C(t) — J duiJiuj) coswt, S(t) = J duJ(u) smut and 
T / — Jo dt\ Jq 1 dt 2 Jq 2 dt 3 , the forth-order coefficients may be written in the following, 



+ [C(t 02 ) sm(u t 03 ) - S(t 02 ) cos(u t 03 )]S(t n ) sm(u t 12 ) 
+ [S(t 3) sm(u t 02 ) - 3C(t 3) cos(oj t 02 )]C (t 12 ) sin(w ii3) 
+ [C(t 03 ) sin(w t 02 ) - S(t 03 ) cos(uj t 02 )]S(t 12 ) sm(u ti 3 ) 
+ \-S(t 03 ) sm(u t 01 ) - C(t 03 ) cos(uj t 01 )]C(t 12 ) sin(w t 23 ) 
+ [-C(t 03 ) sin(cj t i) + S(t 03 ) cos(uJotoi)}S(t 12 ) sin(u t 23 )}, 

S™(t) = 21 J {[S(t 02 ) sm(u t 03 ) + C(t 02 ) cos(u t 03 )}C(t 13 ) sm(u t 12 ) (A5) 
+ [C(t Q2 ) sm(u t 03 ) - S(t Q2 ) cos(u t 03 )]S(t 13 ) sin(u t 12 ) 
+ [S(t 03 ) sm(u to 2 ) + C(t 03 ) cos(u t 02 )]C(t 12 ) sm(u t 13 ) 
+ [C(t 03 ) sm(u t 02 ) - S(t 03 ) cos(u) t 02 )}S(t 12 ) sm(u ti 3 ) 
+ [-S(t 03 ) sin(w ioi) - C(t 03 ) cos(uj toi)}C(t 12 ) sm(u t 23 ) 
+ [-C(*03) sm(u t 01 ) + S(t 03 ) cos(u t 01 )]S(t 12 ) sin^o^)}, 



rf (t) = -8T J {[C(t 13 ) sin(w to3) ± S(t 13 ) cos(u to3)]C(t 02 ) sm(co t 12 ) (A6) 
+ [C(t 12 ) sm(u t 02 ) ± S(t 12 ) cos(u t 02 )]C(t 03 ) sm(u t 13 ) 
=F [S(t 03 )C(t 12 ) + C(t 03 )S(t 12 )] sin(cu t23) cos(^oi)}, 




F (t) = 16T J {[C(t 02 )C(t 13 ) + S(t 02 )S(t 13 )\ sin(a; t 3) sin(u;oti2) (A7) 
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+ [C(t 03 )C(t 12 ) + S(t 03 )S(t 12 )} sin(woto2) sin^ia) 
+ [C{t 0Z )C{t l2 ) - S(t 03 )S(t 12 )} sin(uj t 01 ) sm(u t 23 )}, 



a IV (t) = -8T / {S(t + t 2 )S(t 13 ) sin co (t + t 3 ) sm(to t 12 ) (A8) 
+ S(t + t 3 )S(ti 2 ) smLu (t + h) sm(u ti 3 ) 
+ [C(t 03 )C(t 12 ) - S(t 03 )S(t 12 )} smu (t + ti) sin(o;ot23)} ) 



P IV {t) = 8T / {S(t 02 )S(t 13 ) cosu (t + t 3 ) sm(u t 12 ) 



(A9) 



+ S(t 03 )S(ti 2 ) cosu (t + t 2 ) sin(w ^i3) 

+ [C(t 03 )C(t 12 ) - S(t 03 )S(t 12 )] cosu {t + h) sm(u t 23 )}. 



Appendix B: Derivation of Bloch equation 

According to the definition of Bloch vector bj(t) = Tr[p(t)<jj], we have bj(t) = Tr[p(t)aj]. 
By inserting master equation (2) into it and after some deduction, one can obtain the 
required Bloch equation. For example, for the component equation concerning b x we have, 



b x (t) = -iTr{[H LS (t),p(t)]a x } + Tx{D[p(t)]a x } + 1x{D'[p(t)]a x }. 



(Bl) 



By use of the circulation property of trace operation and the Pauli algorithm, one easily 
get -iTr{[H LS (t),p(t)]a x } = (SL - S+)b y , Tr{D[p(t)]a x } = -±(r_ + T+ + F )b x , and 
Tx{D'[p(t)](T x } = ab x — f3b y . Summing up them, we thus obtain eq.(21). 

The Bloch eqs.(21)-(23) can also be written as the compact vector form, b = Mb + v, 
with the damping matrix M and drift matrix v given respectively by 



M = 



/-i(r_ + r + + r -2a) 





S_ - S + - 13 





-|(r_ + r + + r + 2«) o 

o _ (r _ + r + )y 



(B2) 
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and v T = ^ 0, 0, T + — T_ ^ • Note that the damping matrix is in block diagonal form. 
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FIG. 1. Evolution of the time-dependent coefficients. The dot-dash line, dot line and solid line 
in (a),(b) and (c) correspond to respectively the evolutions of T_, r + and To, while the solid and 
dot lines in (d) refer to the evolutions of a and (3. Where we choose ujq = IOO70, and A = O.270, 
A = 270 for (a); A = 570, A = 5O70 for (b); A = 4OO70, A = IO70 for (c). The parameters for (d) 
are the same as that of (a). 
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FIG. 2. Time evolution of a(t), with the solid and dot lines corresponding to respectively eqs.(25) 
and (39). The parameters in (a),(b) and (c) are set to be in accordance with that in Fig.l. 
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